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Abstract: The Bose-Einstein condensation of a dilute gas of rubidium-87 atoms was achieved by 
cooling a small number of atoms in a magnetic trap. The effective potential of the trap is to lowest order 
harmonic and under these conditions we estimate the critical temperature and the heat capacity of an ideal 
gas of Bose particles. Previously this has been estimated within the semiclassical approximation. However, 
' because of the small number of atoms, the resulting T c deviates slightly from our result obtained by exact 

numerical calculations. For the gap of the heat capacity at the transition point the deviation is 20%. 
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In the remarkable experiment |1|] the rubidium vapor is magnetically trapped and evap- 
oratively cooled below 170 nK where a finite fraction of the particles is observed to occupy 
the ground state. When further lowering the temperature this fraction increases abruptly, 
signaling a Bose-Einstein condensation. The gas is very dilute and hence well approximated 
by an ideal gas. The confining magnetic field is a time-averaged, orbiting potential whose 
effective potential to first order is a three-dimensional anisotropic harmonic oscillator. The 
Bose-Einstein condensation of an ideal gas in such a potential has been treated earlier in 
the semiclassical approximation j|. We will first shortly review a part of this work using a 
slightly different approach leading to the same results. 

The number of particles of an ideal Bose-gas is given by 

n = n ° + T. j„J,*t_ 1 (i) 

where N is the number of particles occupying the ground-state, \i is the chemical potential 
and the sum runs over all one-particle states. The semiclassical approximation consists of 
approximating this sum by the phase-space integral 
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In our case the classical Hamiltonian is given by H = p 2 /2m + \muj 2 [x 2 + y 2 ) + \muj 2 and 
the integral can be evaluated at the transition point where \i = Nq = 0. This gives the 
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critical temperature 

where £(3) = 1.202. Similarly the internal energy is given by 

, v f ci 3 x<i 3 p # 

( } = 7 (2tt^) 3 e( H -^)/ fcr - 1 ( } 

The heat capacity C = dE/dT is seen to be a function of dfi/dT. The latter can be 
determined by taking the derivative of Eq.(0) with respect to the temperature 0. This 
enables us to calculate the heat capacity below and above T c : 
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C(T-) = — ——Nk = 10.81Nk (5) 
V c ' 15C(3) V ; 

AC = C(T C 7) - C(T+) = ^^lNk = Q.57Nk (6) 

The heat capacity is seen to exhibit a gap at the transition point, unlike the ideal gas confined 
in a box. 

In the experiment the number of particles was N = 20000 and the frequencies of the 
oscillators where given by u z /2tt = 120Hz and cu = u z /y/8. Eq.fl3|) then gives the transition 
temperature T c = 73.5 nK. However, in a recent preprint the frequency is corrected to 
uj z /2tt = 208Hz, and will use this value in our numerical estimates. This frequency gives 
a transition temperature T c = 127.4 nK which is closer, but still not close to the observed 
transition temperature at about 170 nK. 

To improve this method one may include interactions using the mean-field approach, 
leading to corrections dependent of the scattering length aH, which for the rubidium gas 
is determined to be in the range 85ao < a < 140a , where ao is the Bohr radius [|]. These 
corrections lowers the transition temperature to Tc = 125.6 for a = llOao- 

Another improvement is to go beyond semiclassical approximation which is valid just 
for kT 3> Tiu. In order to fulfill this requirement at T c the number of particles must be 
quite large since kT c /hu ~ iV 1 / 3 . However, when %uj is comparable to kT in size, the exact 
sum of Eq.(|l|) converges quite fast, and it is possible to evaluate the sum over the quantized 
oscillator-states numerically. Estimating the critical temperature in this way gives a 3% 
correction and T c = 123.4 nK. This correction is seen to be twice as large as the one due to 
interactions. 

After determining T c the internal energy and specific heat may be evaluated by replacing 
the phase-space integrals leading to Eqs.(|5f|D by a finite number of terms of the exact sums. 
This gives 

C(T~) = 10.20iV£; (7) 
AC = 5.28iVfc (8) 
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the latter differing 20% from the semiclassical result, mainly due to the shifted transition 
temperature. 

The semiclassical approach thus turns out to be inaccurate when dealing with such a 
small number of particles and instead the exact quantum sums should be evaluated when 
calculating thermodynamic quantities. The estimated transition temperature deviates from 
the observed. Probably the dilute rubidium gas is well approximated by an ideal gas, at least 
before the condensation takes place. So the crucial point is how well the harmonic oscillator 
approximates the confining potential and if so, how accurate the oscillator strengths may be 
determined. 
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